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Abstract 

Camanho, Edelstein, Maldacena and Zhiboedov have shown that gravitons can experience a negative 
Shapiro time delay, i.e. a time advance, in Einstein-Gauss-Bonnet theory. They studied gravitons 
propagating in singular “shock-wave” geometries. We study this effect for gravitons propagating in 
smooth black hole spacetimes. For a small enough black hole, we find that gravitons of appropriate 
polarisation, and small impact parameter, can experience time advance. Such gravitons can also exhibit 
a deflection angle less than tt, characteristic of a repulsive short-distance gravitational interaction. We 
discuss problems with the suggestion that the time advance can be used to build a “time machine”. In 
particular, we argue that a small black hole cannot be boosted to a speed arbitrarily close to the speed 
of light, as would be required in such a construction. 


1 Introduction 

Lovelock theories of gravity [T] are an interesting alternative to General Relativity in d > 4 dimensions. It 
is interesting to ask how classical properties of these theories differ from GR, in particular whether these 
theories are pathological in some way. 

In GR, causality is determined by the null cone of the metric. It has been known for a long time 
that this is not the case in Lovelock theories. In general, causal properties of a PDE are determined by its 
characteristic hypersurfaces [2]. For example: sharp features of solutions (e.g. a discontinuity in curvature or 
in the 100th derivative of curvature) propagate along such surfaces; for high-frequency waves propagating in 
a background spacetime, the surfaces of constant phase are characteristic hypersurfaces of the background; 
given some initial data surface E, and a region 17 C E, the region of spacetime where the solution is 
determined by data specified on 17 is bounded by an ingoing characteristic surface emanating from 517. 

In GR, a hypersurface is characteristic if, and only if, it is null w.r.t. the metric. However, in Lovelock 
theories, the location of characteristic hypersurfaces depends on both the metric and the Riemann tensor. 
Generically, such surfaces are non-null w.r.t. the metric: gravity can propagate faster, or slower, than 
light Elll]. Furthermore, different graviton polarizations generically propagate along different characteristic 
hypersurfaces (i.e. at different “speeds”) so such theories are “multirefringent”. 

Another important difference between GR and Lovelock theories has been demonstrated recently [^. 
In GR, the equation of motion is always hyperbolic: it has the character of a nonlinear wave equation 
propagating d{d — 3)/2 degrees of freedom in d spacetime dimensions. For small curvature, the same is 
expected to be true in Lovelock theories (although this has not been proved). However, the equation of 
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motion can be non-hyperbolic when the curvature is large. Under such circumstances, the initial value 
problem is no longer well-posed. Consider initial data, consisting (as in GR) of a (d — l)-manifold E, a 
Riemannian metric hab on E and a symmetric tensor Kab (the extrinsic curvature) on E. From this data 
one can determine whether or not the equation of motion is hyperbolic on E. If it is not hyperbolic then, 
generically, there will exist no solution of the equation of motion arising from such initial data: one cannot 
“evolve forward in time”. Solutions may exist in non-generic situations e.g. for analytic initial data the 
Cauchy-Kowalevskaya theorem may guarantee local existence of a solution. However, this is infinitely fine- 
tuned: if one smoothly perturbs the initial data within a compact region then generically there will exist no 
solution of the equation of motion arising from the perturbed data. 

As an example, consider Einstein-Gauss-Bonnet (EGB) theory (a Lovelock theory). This theory admits 
static, spherically symmetric black hole solutions [B]. Since EGB theory has a dimensionful coupling constant, 
we can talk about “small” and “large” black holes. Ref. [S] showed (using results of Refs. 13 [S]) that, for 
a small black hole with d = 5 or 6, the equation of motion is non-hyperbolic in a region outside the event 
horizon. Hence such a black hole solution must be rejected as unphysical: it is not the Gauchy evolution of 
initial data for which the equation of motion is hyperbolic Q 

As we explain in the Appendix, a hyperbolic PDE defines a natural division of vectors into “timelike”, 
“spacelike” and “null” such that causal propagation occurs along non-spacelike directions. In GR, these 
definitions coincide with the notions of “timelike” etc defined using the metric but this is not the case in 
Lovelock theories. Since the causal structure defined by the equation of motion is more important than that 
defined by the metric, we will adopt the following convention: 

Throughout this paper, the words “timelike”, “spaeelike”, “null” or “causal” refer to the causal structure 
defined by the equation of motion, not to the causal structure defined by the metric. If we wish to refer to 
the latter then we will say “timelike w.r.t. the metric” etc. 

There is another restriction on initial data which needs to be mentioned. Given a hyperbolic PDE, for 
a well-posed Gauchy problem, initial data must be chosen so that the initial data surface E is spacelike (in 
the sense just defined). In GR, this is the same as being spacelike w.r.t. the metric but this is not the case 
in Lovelock theories. One can find examples of initial data for which E is everywhere spacelike w.r.t. the 
metric but not everywhere spacelike. Evolution of such initial data is not a well-posed problem. We will 
discuss this further below. 

In summary, there are several restrictions on legitimate initial data in Lovelock theories: data must be 
chosen so that (i) the constraint equations are satisfied; (ii) the equation of motion is hyperbolic; and (iii) 
the initial data surface is spacelike. 

The present paper is motivated by recent work of Gamanho, Edelstein, Maldacena and Zhiboedov (GEMZ) 
[9]. They investigated the Shapiro time delay for gravitons in EGB theory. In this theory, linearised gravity 
propagates at the speed of light in a flat background but can travel faster than light in a curved background. 
Hence one might expect that it is possible for a gravitational signal to travel through a curved spacetime 
“faster” than it would in flat spacetime. GEMZ showed that this is indeed the case. They considered the 
Aichelburg-Sexl (AS) “shock-wave” spacetime |l()l ITT] . This is an exact solution of EGB theory which GEMZ 
interpreted as describing the gravitational field of a high-energy particle e.g. a graviton. They showed that a 
second (“test”) graviton that scatters from the first graviton at small impact parameter can indeed experience 
a negative Shapiro time delay, i.e., a time advance. They also showed that the same result can be obtained 
from a scattering amplitude calculation. 

The AS solution is singular: it has delta function curvature localised on a null hypersurface, with the 
amplitude of the delta function diverging along a line within this hypersurface (the worldline of the high- 
energy graviton). Given that the solution is singular, one must ask whether it is physical. The usual reason 
why the AS solution is regarded as physical is because it can be obtained as a limit of smooth black hole 
solutions. The limit involves boosting a black hole and taking the boost to infinity whilst scaling the mass 

^Of course, one could regard such a black hole solution as the evolution of initial data for which the equation of motion is 
not everywhere hyperbolic e.g. a surface of constant t. However, as just explained, the violation of hyperbolicity implies that 
this is infinitely fine-tuned; if one perturbs such initial data then, generically, it will be impossible to evolve the perturbed data. 
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to zero, keeping the energy fixed. Hence one can “regulate” a singular AS solution by replacing it with a 
small, highly boosted, black hole. 

This motivation for regarding the AS solution as physical is fine in GR but it suffers problems in EGB 
theory. First, for d = 5, 6, small black hole solutions are unphysical, as explained abovej^ There is also a 
problem for d > 7. Gonsider initial data describing a boosted black hole. We will show that, for a small 
black hole, if the boost is too large then such initial data is not everywhere spacelike (in the sense defined 
above). This implies that if one makes a generic smooth perturbation to this initial data then there will not 
exist a solution to the equation of motion arising from the perturbed data. Therefore such initial data is 
unphysical: it cannot arise in Gauchy evolution without infinite fine-tuning. Hence, in EGB theory, a small 
black hole cannot be boosted arbitrarily close to the speed of light: there is a speed limit. This suggests that 
the AS solution may not be physical in EGB theory. Therefore it seems desirable to have an independent 
confirmation that time advance is possible. 

We will study the Shapiro time delay for gravitons in static, spherically symmetric, black hole solutions 
of EGB theory. Another motivation for doing this is that, in general, there is no gauge-invariant definition 
of the Shapiro time delay m- The problem is that there is no gauge-invariant way of identifying points of a 
curved spacetime with those of Minkowski spacetime. However, this problem can be overcome for spacetimes 
with suitable symmetries. In particular, for static, spherically symmetric spacetimes, it is possible to define 
the Shapiro time delay unambiguously iniB The idea is to consider a spherical cavity surrounding the black 
hole and calculate the (proper) time it takes for a graviton to cross the cavity, and then compare with the 
time it takes to travel between the corresponding points of a spherical cavity in fiat spacetime. 

Our main result is to confirm that a time advance for gravitons is indeed possible for small black hole 
solutions of EGB theory. This occurs for gravitons encountering the black hole with appropriate polarisation, 
and an impact parameter comparable to the length scale set by the GB coupling, just as GEMZ found for the 
AS spacetime. Furthermore, we find that such gravitons can undergo a deflection through an angle less than 
TT. Both of these features indicate that certain graviton polarizations experience a repulsive gravitational 
interaction at distances comparable to the scale set by the coupling constant, which is confirmed by examining 
the “effective potential” which determines graviton trajectories. Glose to the black hole horizon, gravity 
becomes attractive again, and the deflection angle is greater than tt for very small impact parameter. 

Time advance can occur when the cavity is not too large. For a very large cavity, the time delay is positive. 
This is because the graviton suffers a deflection, which makes its path longer than the corresponding (straight 
line) path in flat spacetimej^ This results in a contribution to the time delay which grows in proportion to 
the cavity radius, ensuring that the delay is positive for large enough cavity radius. However, there are two 
special values for the impact parameter for which the deflection angle is exactly tt, i.e., the graviton suffers 
no net deflection. For these, the time delay is finite in the limit of infinite cavity radius. For the larger of 
these values, the time delay is positive. However, for the smaller value it can be negative. So this is an 
example of time advance for a graviton incident from, and returning to, infinity. 

GEMZ argued that the possibility of a time advance implies that EGB is unviable as a classical theory. 
They employed an argument which asserts that when one has superluminal propagation in a Lorentz covariant 
theory one can construct a time machine, i.e., a spacetime with closed causal curves (where “causal” is 
defined as above). Arguments of this sort have been applied to various non-gravitational field theories in 
Ref. [14]. However, this kind of argument has been criticised by Geroch [15], who notes that the existence of 
a causally pathological solution is not enough to reject the theory (see also Ref. (Hj). After all, GR admits 
causally pathological solutions e.g. Minkowski spacetime with a periodic time direction. The existence of 
such solutions is not regarded as problematic because one cannot “form” them, i.e., they do not arise as 
the Gauchy evolution of initial data. We will argue that the “time machine” constructions of Ref. m and 
GEMZ also do not arise as the Gauchy evolution of legitimate initial data. The reason is that the initial 
data surface is not everywhere spacelike (in the sense defined above). Hence it is not a well-posed problem 
to specify a solution in terms of such data: either no solution exists, or it is infinitely fine-tuned. 

^In fact, for d = 5, there is a mass gap so a black hole cannot have arbitrarily small mass. 

®The large symmetry of the AS solution implies that the time delay should be unambiguous in this case too. 

^CEMZ evaded this effect by considering a graviton propagating between two AS shocks, resulting in zero deflection. 
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If the “time machine” construction doesn’t work then is there some other reason to regard EGB theory, or 
more general Lovelock theories, as “badly behaved” when compared to GR? There is no proof that the initial 
value problem for Lovelock theories is well-posed, even for initial data for which the equation of motion is 
hyperbolic and the initial surface is spacelike e.g. small curvature initial dataj^ If the initial value problem 
is not well-posed then these theories don’t make sense classically. Assuming that the initial value problem is 
well-posed. Ref. uni argued that Lovelock theories, unlike GR, will suffer from shock formation similar to a 
compressible perfect fluid. Since different parts of a wavepacket can travel with different speeds, one could 
construct a smooth wavepacket for which the back travels faster than the front, leading to wave steepening 
and, eventually, a shock. The curvature would diverge at a shock. Ref. [19] argued that this would be a 
naked singularity. So this is an example of a pathology that afflicts Lovelock theories but not GR. However, 
it was argued to be a “large data” effect, i.e., it may not occur for small curvature initial data. Indeed, 
Ref. [T2| argued that Minkowski spacetime is likely to be nonlinearly stable in Lovelock theories. 

This paper is organised as follows. We first give an informal review of the notions of hyperbolicity and 
causality required in EGB theory (Section]^. In Section we review the static, spherically symmetric, 
black hole solutions of EGB and derive a “speed limit” for small black holes. In Section]^ we will investigate 
the Shapiro time delay and deflection of gravitons propagating in a static, spherically symmetric, black 
holes solution of EGB. In Section]^ we discuss the “time machine” constructions of Ref. [IT] and GEMZ. 
Appendix gives a more technical review of how a hyperbolic PDE defines a notion of causal structure. 
Appendix]^ gives details of our perturbative calculation of the time delay and deflection angle. 

2 Hyperbolicity and causal structure 

In this section we will give a brief review of the definitions of hyperbolicity and causality that we will need. 
A more technical review is presented in Appendix Hyperbolic PDEs describe phenomena that exhibit 
finite propagation speed. A hyperbolic PDE defines a causal structure: given a spacetime point p, one can 
define the causal past of p as the set of points q for which a change in the solution at q would result in a 
change in the solution at p. Similarly the causal future of p is the set of points q where the solution at q 
which depend on the solution at p. In general, the causal past/future depend on the solution itself. 

In general, the causal structure is determined by the coefficients of the highest order derivatives in 
the PDE. The causal structure of a PDE is closely related to its characteristic hypersurfaces. Loosely 
speaking these are hypersurfaces on which all of the highest order derivatives cannot be determined in terms 
of all the lower derivatives by the equations of motion. If one has a solution which is discontinuous in the 
highest derivatives (but continuous in the others) across a hypersurface, then that hypersurface is necessarily 
characteristic. Hence, discontinuities in the solution propagate along characteristics. More “physically”, if one 
considers high frequency waves propagating on top of a “background” solution then the surfaces of constant 
phase are characteristic hypersurfaces of the background solution. 

There is an algebraic criterion to determine whether or not a hypersurface is characteristic. Given a point 
p and a covector ^ at p one can define the characteristic polynomial Q{p,f), a homogeneous polynomial in ^ 
with coefficients determined by the coefficients of the highest derivative terms in the PDE. A hypersurface 
S is characteristic if, and only if, its normal one-form f satisfies: 

Q(p,O = 0 (2.1) 

for any p G E. At each point, the set of one-forms satisfying this equation defines a cone in the cotangent 
space, known as normal cone. 

In GR, the normal cone is simply the light cone in the cotangent space. In Lovelock theories, the 
normal cone is generically non-null w.r.t. the metric mm- In GR, the d{d — 3)/2 degrees of freedom 
of the metric all propagate with the same speed. However, in Lovelock theories these degrees of freedom 
generically propagate at different speeds so the normal cone is expected to be composed, generically, of 

®See Ref. 11711181 for discussion of the first steps in constructing such a proof. 
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d{d — 3)/2 distinct “sheets”]^ For linearised perturbations of a solution, this corresponds to the statements 
that different graviton polarizations propagate with different speeds: Lovelock theories are multirefringent. 

The definition of hyperbolicity is reviewed in the Appendix. Roughly speaking, this amounts to requiring 
the “correct number” of propagating degrees of freedom at each point. For the equation of motion to be 
hyperbolic, we require that the normal cone has d{d — 3)/2 sheets. Obviously this isn’t true in GR, so 
one has to allow for degeneracy i.e. the possibility that multiple graviton polarizations propagate along 
the same sheet. We also require that the interiors of these sheets have a non-trivial intersection. We then 
define a covector to be timelike if it lies in the interior of this intersection, causal if it lies in the closure of 
this intersection (and is non-zero), and spacelike if it lies outside this closure. A hypersurface is said to be 
spacelike if its normal covector is timelike everywhere. Hence if the theory is hyperbolic then (i) we have 
the right number of degrees of freedom and (ii) we can find surfaces that are “spacelike for all degrees of 
freedom”. 

The dual of the normal cone, is another cone, this time in the tangent space, known as the ray cone. 
In the case of the linear wave equation or General Relativity this is simply the null cone defined by the 
spacetime metric. In a Lovelock theory it will generically be composed of d{d — 3)/2 distinct sheets. The 
ray cone is generated by hicharacteristic curves. These are determined, in local coordinates x^, by solving 
the Hamiltonian equations 

ds ds dx^ 

with ^(0) chose so that (5(x(0),^(0)) = 0. For the wave equation, or for GR, bicharacteristic curves corre¬ 
spond to null geodesics. For high frequency waves, one can use the geometric optics approximation to show 
that high-frequency gravitons follow null geodesics. Similarly, in a Lovelock theory, high-frequency graviton 
worldlines are bicharacteristic curves (see e.g. [TS]) so we will calculate the Shapiro time delay for such 
gravitons by determining these curves. 

We define the causal cone as the convex hull of the sheets defining the ray cone. We will then define a 
vector to be causal if it lies in its closure and spacelike if it lies in its exterior. Note that a hypersurface is 
spacelike iff all tangent vectors are spacelike. 

To illustrate the above definitions, consider the following example. For some non-generic spacetimes 
(possessing appropriate symmetries) the characteristic polynomial factorizes into a product of quadratic 
factors: 

Q{p,o = {Gi{pr'‘UbYKG2{pr'‘UbY^ ■ ■ ■ ( 2 . 3 ) 

For example, this happens for any Ricci fiat spacetime with a type N Weyl tensor, and also for the static 
spherically symmetric spacetimes that we will consider below . Each of the Gf’ can be interpreted as an 
(inverse) effective metric. A hypersurface is characteristic iff its normal is null w.r.t. one of these effective 
metrics. The normal cone is simply the union of the null cones of Gf^. In this case, the theory is hyperbolic 
iff (i) there are the “right number” of effective metrics, i.e., one for each degree of freedom (allowing for 
the possibility that one effective metric may correspond to several degrees of freedom); (ii) each Gf' has 
Lorentzian signature; and (iii) the interiors of the null cones of Gf have a non-trivial intersection. This is 
the case for any Ricci fiat type N spacetime [S]. However, as we will see below, for certain small black hole 
solutions, it turns out that one of the Gf can fail to be everywhere Lorentzian, so the theory is not hyperbolic 
in such spacetimes [5]. For both examples, when the theory is hyperbolic, it turns out that the null cones 
of Gf form a nested set, so causality (for covectors) is determined by the effective metric corresponding to 
the innermost null cone. We define Gjat to be the inverse of Gf. Then the ray cone is the union of the null 
cones of Gjab and causality (for vectors) is determined by the effective metric with the outermost null cone. 
The bicharacteristic curves are simply the null geodesics of the G/ab- 

®This has been proved only for the class of Ricci flat type N spacetimes [^. 
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3 Einstein-Gauss-Bonnet black holes 

3.1 Spherically symmetric EGB black holes 

The equation of motion of Einstein-Gauss-Bonnet (EGB) theory is obtained by varying the action 


^ ^ [R + Xgb{R^ - (3.1) 

where Xgb is the GB coupling constant and we have set G = 1. This theory admits static, spherically 
symmetric solutions with metric |3]: 

9 = + /(r)"Mr^ + r'^gsd-2 (3.2) 

where gsd -2 is the standard round metric on the unit {d — 2)-sphere and 

/(r) = 1+—(l-q(r)) (3.3) 

a 

«(’■) = \/i + ^ 

a = 2{d - 3){d - 4)Xgb (3.5) 


and we assume a > 0 since for a < 0 the metric has a naked singularity The parameter g is proportional 
to the ADM mass M: 

_ (d - 2 )Vo1(S^-2) 

M=^-- ’-g. 3.6 

IGtt 

The event horizon is a hypersurface r = r// where th is the largest root of 

9 = (xr + —'j ■ (3.7) 

The coupling a has dimensions of length squared so EGB theory has a length scale ^/a. We will say that a 
black hole is “small” if th ^ y/a and “large” if th ^ y/a. Equivalently, for d > 5, a black hole is small if 

g < (3.8) 

and large if g^ Note that black holes with arbitrarily small mass do not exist for d = 5 because 

there is a mass gap (however tr can be arbitrarily small): g > a/2. The function q varies over a length 
scale 

L = {ga)x^. (3.9) 

For r ^ L we have 

(3-10) 

i.e., the solution reduces to the higher-dimensional Schwarzschild solution. For a large black hole rn ^ L 
so this approximation is valid everywhere outside the horizon. 

priori, the sign in front of q in the expression for / is arbitrary. We choose the negative branch, corresponding to 
asymptotically flat solutions. 

®Since some modes of gravitational perturbations can travel faster than light, it could be possible, a priori, that such 
perturbations could escape the black hole region (as defined by the causal structure given by the physical metric). However it 
was shown in [20||5] that a Killing horizon is always a characteristic hypersurface for all graviton polarizations, excluding this 
possibility when the event horizon is a Killing horizon, as is the case here. 
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3.2 Effective metrics and bicharacteristic curves 


Characteristic hypersurfaces of the above solution were determined in Ref. |5] . The symmetries of the solution 
imply that the characteristic polynomial factorizes into a product of quadratic factors, each associated to an 
“effective metric” as discussed above. A hypersurface is characteristic iff it is null w.r.t. one of these effective 
metrics. The explicit form of the effective metrics was determined by considering linear perturbations of 
such solutions. Such perturbations can be classified into scalar (S), vector (V) and tensor (T) types w.r.t. 
the spherical symmetry. For each type, one can obtain a “master equation” mi and from these one can 
read off the effective metric for that type. Hence the characteristic polynomial factorizes as 

Q(p, 0 = {GfUb^ {G’ifUdY- {gG Ufr (3.11) 


where Ps,Pv,Pt denote the number of degrees of freedom of each type of modes. Viewing etc as inverse 
metrics, the corresponding metrics are given by 


GA = -f{r)dY +fir)-^dY + 


for certain smooth functions caY) given by [5] 

cs{r) =3(1- Ar) + ( 1 - 


cv{r) = A{r), 

CtY) = - (1 + ) A{r) - ( 1 - 


d-4 



^2 


CAir)' 

3 \ 

1 

d-2) 

Air) 

1 ^ 

\ 1 

d — 4y 

1 Air) 


-3 1- 


d-2 r 


+ 3. 


where 


Air) = q{r)-^ U + :^ + U " 


1 1 


2 d-3j \2 d-3, 

It is convenient to take the index A G {0, S, V, T} where 0 refers to the physical metric, i.e.. 


(3.12) 


(3.13) 

(3.14) 

(3.15) 

(3.16) 


Coir) = 1 . 


(3.17) 


For a large black hole, the functions c^(r) are positive everywhere outside the horizon. This ensures that the 
effective metrics have Lorentzian signature, and their null cones form a nested set, with the outermost cone 
(in the tangent space) corresponding to the effective metric with the largest value of caY). (The physical 
metric Go is included in this nested set.) This ensures the hyperbolicity of the theory, and causality is 
determined by this outermost null cone. If ca > 1 then the associated modes propagate faster than light [5]. 

For d = 5,6, for small enough r^j, it turns out that one of the ca('c) vanishes at some value r = > r[{ 

and becomes negative for r < [^. The corresponding inverse effective metric G'Y is smooth at r = r* 

but becomes degenerate there. For r < r* it has Lorentzian signature, but with the opposite overall sign. 
This implies that the theory is non-hyperbolic in such black hole spacetimes. Therefore small black holes 
are unphysical for d = 5, 6. This does not occur for EGB theory with d > 70 


Note that the function ca are determined entirely by the length scale L defined by (3.91. If r ^ L then 


caY) = 1 + 2/3a ( — 

' r 


d-l 



where the constants Pa given by 

(d-i) 


/3s = - 


id-sy 


yv = -T; 


lid-1) 


2 id-3)' 


Pt = 


2(d-l)N 


(d-l) 


(d-3)(d-4)’ 


/3o = 0. 


(3.18) 


(3.19) 


^However, it does occur for more general Lovelock theories with d > 7. Generically, it occurs when the equation of motion 
includes the highest order Lovelock term [5]. 
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We see that < 1 at large r so the scalar and vector polarizations propagate slower than light in this 

region. However ct > 1 at large r so tensor polarizations propagate faster than light at large r. Hence 
causality at large r is determined by the effective metric for the tensor modes. 

We will now prove that cg < 1 and cy < 1 everywhere. Since q(r)~‘^ is monotonically increasing we 
see that also A{r) is monotonically increasing. We also have ^(oo) = 1. It follows that .4(r) < 1 hence 
cy(r) < 1. Now we look at cs{r). We have: 


c'sir) = 


3 1 - 


d-2 


- 1 - 


d-2 


Air) 


-2 


A'{r). 


(3.20) 


Since .4(r) is monotonically increasing, the sign is determined by the terms in parentheses. For d = 5 this 
is constant and positive, hence > 0 so csir) < 05 ( 00 ) = 1. For c? 7 ^ 5, the expression in parentheses is 
negative at small r and positive at large r. Hence, starting from r = 0, csir) decreases to a minimum and 
then increases monotonically with r. Hence csir) < max{c 5 (oo), 05 ( 0 )} = 05 ( 00 ) = 1. 

The same argument allows us to determine an upper bound for ct- We have: 


c'rir) 


1 + 


d-4 


- 1 - 


1 


d-4 


Air 


,-2 


A'ir). 


(3.21) 


If d = 5 then the expression in square brackets is constant and positive so we see that ct is monotonically 
decreasing hence 

CT(r) < ct(0) = 3 (d=5). (3.22) 

If d > 5 then ct has a maximum at vq where 


A(ro) = 


'1- l/(d-4) 
l + l/(d-4) 


and hence 


crir) < criro) = 3 - 2 Jl - 


(d-4)2 


Note that the RHS is greater than 1 . 


(3.23) 


(3.24) 


3.3 Speed limit for small black holes 

We will now consider the effect of boosting one of these black holes. To construct initial data describing a 
boosted black hole, we can consider the data induced on a boosted hypersurface in the black hole spacetime. 
Such a hypersurface is spacelike w.r.t. the metric for any boost velocity v such that |v| < 1. However, since 
the null cone for the tensor modes can lie outside the null cone of the physical metric, it is possible that, for 
|u| close to 1, the hypersurface may fail to be everywhere spacelike w.r.t. the tensor effective metric. This 
implies that it will fail to be spacelike in the sense defined in the Introduction and hence it would not be a 
valid initial data surface. We will now show that this is indeed what happens. 

First we introduce an “isotropic” radial coordinate f defined by 



d log r 1 

dr r\/J 

(3.25) 

so that the physical metric is 

g = -fdf + H idf^ + f‘^gsd- 2 ) 

(3.26) 

where 

2 




(3.27) 
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For r ^ L we can use the approximation (3.10) to obtain 


1 - 


2(d- 3)H-3 


(3.28) 


and hence 


/ 


1 - 


^d—3 


H^l + 


(d-3) 


xd—3 ' 


(3.29) 


To construct initial data describing a boosted black hole we convert to Cartesian coordinates so that 
= f cos 01 etc (where 61 , 62 , ■■ ■ are the angles on and then perform the Lorentz transformation 


—vt'), t = "fit'— vx^ ), 7 = (1 — u^) 


(3.30) 


We now consider the data induced on a surface of constant t'. By inverting the Lorentz transformation, we 
see that this is the same as the data induced on a surface of constant t + vx^, i.e., a surface of constant 
t + uf COS 01 . Let S be such a surface. Define a 1-form ^ normal to E: 


^ = dt + v cos 0 idf — vfs\i\- 6 id 6 i. 


(3.31) 


We want to take the data induced on E as initial data describing a boosted black hole. To do this, we must 
check that E is spacelike in the sense defined in the Introduction. This is equivalent to E being spacelike 
w.r.t. all of the effective metrics, i.e., ^ must be timelike w.r.t. to all of the effective metrics. To investigate 
whether this is this case, consider the norm of ^ w.r.t. Ga at 61 = 7r/2: 

= -r' + (3.32) 

Now assume that the black hole is small and consider the region L ^ r ^/a where . 

Using our expansion for ca we then obtain 

G'X^^^.U9,=./2 - -(1 - v'^) + 2 ^av^ . (3.33) 

For the scalar, vector and physical metrics we have /3 a < 0 so the RHS is always negative. However, for the 
tensor effective metric we have Pt > 0 and hence if v is too close to 1 then the second term above, although 
small, will overwhelm the first and the RHS will be positive, i.e., ^ will be spacelike w.r.t. Gt and hence E 
will not be everywhere spacelike. 

As discussed in the Introduction, it is not a well-posed problem to evolve initial data if E is not spacelike. 
Of course, we know that this particular data on E can be evolved - the resulting solution is just the black 
hole solution described above. However, the lack of well-posedness implies that this procedure is infinitely 
fine-tuned: if we make a generic (smooth) perturbation to the initial data on E (for v very close to 1) then 
it will not be possible to evolve the perturbed data either forwards or backwards in time. Hence there is a 
speed limit for small black holes: they cannot be boosted to velocities arbitrarily close to the speed of light. 

One might criticise this argument on the grounds that there is no unique way to boost a black hole. One 
could consider a different surface which is asymptotic to E but differs in the region L r y/a in which E 
can fail to be spacelike. However, note that in our argument, the physical metric is actually flat to the level 
of approximation used because we neglected terms of order The boost used above is a symmetry of 

this flat metric. Therefore our surface E conforms to the usual idea of a boosted hypersurface in the region 
relevant for the above argument. Of course the effective metrics are not flat to this level of approximation 
because they include larger terms, of order — iL/rY~^. 

It is interesting to determine the critical value of v below which E is spacelike. For E to be spacelike 
w.r.t. the tensor metric we need ^ to be timelike everywhere, which is true ifif^ 

^2 

< ^^max = min (3-34) 

r^JCT 

used 0i = 7r/2 to derive this. It is not hard to show that other values of 6i gives less stringent constraints. 
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where the minimum is taken over, say, r > th- For a small black hole, the minimum is achieved for 
L r a/q") which / « 1 and r ^ r hence 


1 

max ct 


3-2,/I- 


< 1 . 


(d-4)2 


(3.35) 


This is the speed limit for an arbitrarily small black hole. More generally, Umax will depend on the mass of 
the hole, with Umax —t 1 for a large black hole. 

It is natural to ask what would happen if one attempted to accelerate a small black hole to a speed 
greater than Umax- As emphasized in Ref. m, one would have to specify the details of how one would 
attempt to achieve this acceleration using only the fields present in the theory. Perhaps one could set up 
initial data consisting of several black holes in the hope that a “gravitational slingshot” effect could be used 
to accelerate a small black hole to a speed greater than Umax- However, as we will see in more detail below, 
the gravitational interaction associated with small black holes in EGB is very different from GR so there 
is no reason to believe that a small black hole in such a system would behave in the same was as it would 
in GR. Whatever the system does, it will not result in a small black hole moving with a speed arbitrarily 
close to the speed of light at some later time. This is because an “instant of time” corresponds to a spacelike 
hypersurface and the argument above excludes the possibility of a small black hole moving arbitrarily close 
to the speed of light on such a surface. 


3.4 Graviton trajectories 

As discussed above, characteristic hypersurfaces are generated by bicharacteristic curves and, in the present 
case, these are simply the null geodesics of the effective metrics. Hence, in the geometric optics approxima¬ 
tion, the worldlines of high-frequency gravitons are null geodesics of the effective metrics. We will need to 
determine these geodesics in order to calculate the time delay. 

Gonsider a null geodesic of GAab- Introducing polar coordinates (01,02, - - - ,0d-3,</>) on spherical 

symmetry allows us to assume that the geodesic is confined to the equatorial plane 0i = - - - = dds = 7r/2. 
Associated to the Killing fields d/dt and d/dcf) are the conserved quantities 

E = f{r)i J=-^^ (3.36) 

CA[r) 

where a dot denotes a derivative w.r.t. an affine parameter A. E and J are not physical because they depend 
on the choice of affine parameter. However their ratio is independent of this choice: 


b = 


J 

E 


(3.37) 


and this is the impact parameter of the geodesic. The null condition gives 




where the effective potential is given by 


Vlf{r) = 


fir)cAir) 
2r2 ■ 


(3.38) 


(3.39) 


Plots of the effective potential for some different cases are given in Ref. and also in Figure 

The effective potentials exhibit a local maximum corresponding to an unstable graviton orbit analogous 
to the photon sphere in GR. Hence in EGB there is a distinct “graviton sphere” for each graviton polarisation. 
We will refer to these as the “scalar sphere”, “vector sphere” and “tensor sphere”. (In some cases, it turns 
out that there are two local maxima and a local minimum of the effective potential and hence three graviton 


10 








d=5,r„=l 





d=8,rH=l 



Figure 1: Effective potentials for a black hole with r// = 1 in d = 5,6, 7, 8 dimensions. We fix the Ganss- 
Bonnet conpling Agb = 2. The red cnrve corresponds to the effective potential for photons, i.e., nnll geodesics 
of the physical metric (c = 1). Superluminal propagation (ca > 1) corresponds to an effective potential which 
is larger than that for photons. This happens only for tensor polarizations. The violation of hyperbolicity is 
associated with the region in which one of the effective potentials becomes negative. This happens near the 
horizon for small black holes in five and six dimensions. 


spheres.) If Emax denotes the maximum of the effective potential then a graviton incident from large distance 
will be absorbed by the black hole if 


6 " < 


2K 


(A) 


absorption. 


(3.40) 


We will consider only gravitons with larger impact parameter, which are scattered by the black hole. 


For r ^ L, equations (3.101 and (3.18) imply that the effective potentials have the expansion 


= 


A* 


2j»2 2r'^~^ 


-Pa 


afx 


(3.41) 


The first two terms are familiar from GR: the first is a centrifngal barrier and the second is responsible for 
the deflection of light rays and the time delay of photons. The third term arises from the Ganss-Bonnet 
interaction. For the effects of this term to be non-negligible compared to the second term we need r < y/a. 
Since we have assumed r ^ L this requires L <C y/a, which implies (3.8), i.e., the black hole has to be small 
compared to the GB scale for this term to be important]^ Notice that this term is negative for vectors and 
scalars but positive for tensors. Hence, for a small EGB black hole, tensor-polarized gravitons experience a 


^This is not possible for d = 5, 6 because of the failure of hyperbolicity for small black holes with d = 5, 6. 
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new repulsive interaction for L < ^/a. It is this repulsive interaction that allows for the possibility of a 
time advance. 

For a small black hole, the effective potential also simplifies in the region L ^ r ^/a. In this region we 
can approximate / ~ 1 and hence 


V^^Hr) 


CAjr) 

2r2 


(3.42) 


which depends only on the length scale L. For tensor modes, the RHS typically has a maximum at some 
value r ^ L so we deduce that the “tensor sphere” has r ~ L. Since ca 1, we deduce from (3.401 that 
tensor-polarized gravitons which scatter from the black hole (rather than being absorbed) must have b > Lp^ 
Note also that the presence of this maximum implies that the interaction between tensor-polarized gravitons 
and the black hole must become attractive again for r < L. 


4 Time delay and time advance 

4.1 Photon time delay in GR 

As discussed in the Introduction, there is no gauge-invariant definition of the Shapiro time delay applicable 
to a large class of spacetimes HH. However, for a static, spherically symmetric spacetime, there is an 
unambiguous definition of this quantity m- The idea is to compare the time it takes a photon to travel 
between two points of a spherical cavity with the corresponding time in Minkowski spacetime. In more 
detail, one can introduce coordinates (in 4d GR) so that the metric takes the form 

ds^ = -A{r)dt^ + B{r)dr^ + r^dfl^ (4.1) 

with A, B > 0. Now one can consider a photon trajectory which starts and ends at r = i?, with r < R along 
the trajectory. Using the spherical symmetry, we can assume that the motion is confined to the equatorial 
plane 9 = 7 r/ 2 . Assume it starts at t = to, (f) = (j)o and ends at t = to + At, cj) = (j)o + Acj). The coordinate 
time to traverse the cavity is At. The proper time (according to a cavity observer) is 

At = yjA{R)At. (4.2) 

One can compare this with corresponding quantity in Minkowski spacetime where the trajectory is simply a 
straight line traversing the cavity, which takes proper times 

ArMink = 2 i?sin • (4-3) 

Hence the time delay can be defined as 

D = At - ArMink- (4.4) 

A simple argument (based on [IS]) shows that, in GR, it is impossible to have D < 0 for a large class of 
spacetimes of the above form. Gonsider spacetimes for which A'{r) > 0 and B{r) > 1. Examples of such 
spacetimes include the (positive mass) Schwarzschild solution and perfect fluid stars with positive energy 
density and pressure. Now let A G [0,1] be an arbitrary parameter along the photon trajectory (which has 
r(A) < R). Since the photon trajectory is null we have 

At = [ \lA~^Bf‘^ -(- A~^r‘^(j>'^ dX > — [ \/ -I- dX > — ^ sin ( —(4.5) 

Jo ^ ./A(R^ Jo ^ V 2 ; ^ ^ 

where the first inequality follows from B{r) > 1 and A'(r) > 0, so A(r) < A{R) for r < R. The second 
inequality follows from the fact that the distance in Euclidean space is minimized by a straight line. It 
follows immediately that D >00 

^^Hence the absorption cross-section for (high frequency) tensor-polarized gravitons by a small black hole scales as 
rather than 

If we normalize t so that A{r) ^ 1 as r ^ oo then A'(r) > 0 implies A{r) < 1 everywhere so Ar < At. Hence the time 
delay defined using At instead of Ar also will be positive m- 
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4.2 Time delay in EGB 


We can now calculate the Shapiro time delay for gravitons propagating across a spherical cavity in the 
geometry (3.2 1 . The cavity is taken to be the surface r = i?13 Consider a graviton worldline parametrised 
by A S [0,1] that has r < R and starts and ends at r = R with (/)(0) = 0, ^(1) = Acj). From the fact that this 
world line is null w.r.t. the relevant effective metric, the coordinate time t taken for the graviton to traverse 
the cavity is 

f-l 1^2 7.2 , 

7 ^ + 7 — dX. (4.6) 

fCA 


(4.7) 


At = 

Jo 

One can show that / is monotonically increasinso f{r) < f(R) < /(co) = 1. Hence 

1 ' ' 



At > 





^ dX > -^= 

f CA ./m 



We showed above that cg < 1 and cy < 1. This implies that for scalar or vector polarizations (or for photons, 
which have cq = 1 ) we have 


At > 


1 


^/7{R) 


J .2 _|_ j. 2^2 > sin 

y f{R) 


A(j) 


ArMink 

TtW) 


(4.8) 


where, as before, ArMink is the time it takes a photon (or graviton) in Minkowski spacetime to travel across 
the cavity between the same two points. Converting to proper time we therefore have 


= \/f{R) At > ArMink 


(4.9) 


so the time delay is always positive for scalar or vector polarized gravitons. More physically: gravitons 
with these polarizations travel slower than photons so, since photons experiences a positive time delay, these 
gravitons must also experience positive time delay. 

The story is different for tensor modes: as we have shown, C 7 ’(r) can be larger than one, so one cannot rule 
out the possibility of time advance (e.g. in d = 5 we would have crir) G [1,3] and thus Ar > ArMink/v^)- 
In the next subsection we will show that time advance is indeed possible, in agreement with [S]. 


4.3 Time advance in EGB: perturbative results 


We will now show how one can achieve a negative time delay, i.e., a time advance, for gravitons of tensor 
polarizations in the space time of a small black hole. We will calculate the time delay explicitly. For 
completeness, we will also present results for the time delay for vector and scalar graviton polarizations, and 
also for photons. Consider a graviton trajectory, given by a null geodesic of the relevant effective metric. As 
before, we assume that this starts at t = 0, r = R, cj) = 0 and ends at t = At, r = R, (p = Ap with r < R 
along the trajectory. Let Rq be the minimum value of r along the trajectory. As in GR, Rq uniquely labels 
the trajectory. This is related to the impact parameter b as: 


/(i?o)cA(i?o)' 


(4.10) 


We can compute the proper time and the deflection angle from the geodesic equations: 

A4. = 2b£drcAr) 

^^Note that R is the area-radius of the cavity w.r.t. the physical metric, not w.r.t. an effective metric. 
Use q' = —{d — l){q^ — l)/{2rq) to show that f' > 0. 


(4.11) 

(4.12) 
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As before, this includes results for photons (with cq = 1). Recall that Rq must be larger than the radius of 
the phot on/graviton sphere for the physical/effective metrics. Both the time delay and the deflection angle 
will diverge as Rq approaches this value since the corresponding trajectories will orbit the black hole many 
times. 

We will first calculate the above quantities for a graviton tra.iectory that has Rq ^ L and also Rq ^ 
^i/(d-3)^ and hence r ^ L, r ^ along the trajectory!^ Under these assumptions, b and Rq are 

related by 


& =i?n 


1 + 


^ _ ‘2a/A \ 


R\ 


1-3 


Rl J 


(4.13) 


We will assume also that the cavity radius is large: R^ Rq. The perturbative calculation of the deflection 
angle is explained in Appendix^ The result is 


This has a well-defined limit when the cavity radius is taken to infinity at fixed Rq. Using b : 
the result in terms of b gives 


(4.14) 


Rq to write 


= lim A(() « TT -k 

R^oo 0“ 


V^{d-1) 1 - 


2aPA{d-2)\ r{l) 


62 


{d-l)J 2T{‘^)' 


(4.15) 


This is the analogous to the result for the deflection of light by the Schwarzschild solution in GR. Note that 
A/ao > TT a Pa A 0 SO scalar and vector polarised gravitons, and photons, are always deflected towards the 
black hole. However, for tensor polarised gravitons, since Pt > 0) we see that A^oo < when 


b< 


' 2{d-2)a 

(d-3)(d-4)' 


(4.16) 


This is consistent with our previous assumptions if, and only if, (3.8) holds, i.e., iff the black hole is small. 
Hence, for a small black hole, tensor-polarized gravitons with b obeying ( 4.16| ) (and b Rq ^ L from 
our previous assumptions) are deflected away from the black hole. This is precisely because such gravitons 
experience the repulsive short-distance interaction in (3.41) that we discussed above. 


Using the same approximations as above we And that the time it takes a graviton to cross the cavity is 


At = 2R- 




i?: 


d—A 


TT 1 — 


2aPA (d - 4) \ (d - l)(d - 3) T (|) 


Rl {d-3)J (d-4) 2r(^^^) 


(4.17) 


and R^ Rq implies that At = At to this level of approximation. The corresponding time in flat spacetime. 


with deflection angle Ap is given by (4.3), which can be written as 

' n — Ap 


ArMink = 2R cos 


= 2R 


1 / TT — Ap 

~ 2 [ 2 


so plugging in our perturbative result (|4.14| gives 
ArMink = 2i? 




R, 


+ 0 R 


d-i 


(d_i) r(^) 


Rl (d-l) 


2r(‘^; 


R 


d—3 


rd-1- 


r: 


d-l 


(4.18) 


(4.19) 


(4.20) 


^®For a large black hole, i.e., rn ^ these conditions reduce to Rq ^ vh- For a small black hole, i.e., th they 

d — 5 ^ 

reduce to Rq ):§> a ■ 
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Note that second order corrections to this result grow linearly with R. For these to be small compared to 
the terms that we have retained we need 


^ _ _ ^ ( fill ) = tk _ 

Ro n ’ Ro \L J an 


(4.21) 


i.e. the cavity is large {R/Rq ^ 1) but not too large. If these assumptions are not satisfied then most of the 
trajectory is in a region where spacetime is almost flat and a large positive time delay (proportional to R) 
results simply because, in this flat region, there is a shorter (straight line) path available which remains far 
from the black hole. 

Combining the above results gives the time delay as (using Rq « b) 


D = At - ArMink 


fJ- 

bd-i 


2a/3^(d-4)\ (d-1) r(f) 

62 {d-l)J {d-4)2r{^) 


(4.22) 


We see that scalar or vector polarised gravitons, or photons, suffer a positive time delay as expected. However, 
for tensor polarised gravitons, a negative time delay, i.e. a time advance results when 


6 < 


2a 


(d-3)- 


(4.23) 


Since b ^ Rq, this is consistent with Rq ^ L, as assumed above, only for a small black hole. For such a 


black hole, stated in terms of 6, our assumptions in deriving (4.22) are 


6> L, 


R 6‘^-3 

1 < y- < -: 

b n 


d-l 


(4.24) 


Note that this overlaps the region for which Acpao < i-e., the gravitons that experience a time advance are 

also deflected away from the black hole. Both effects arise from the repulsive term in the effective potential 
discussed above. 

As discussed above, we need to impose an upper bound on R/b to see the time advance since for very large 
R there will be a large time delay, proportional to R, which occurs because the trajectory has undergone a 


deflection. However, for the special case of a trajectory which saturates (4.16), we have A^oo = tt, i.e., there 


is no net deflection (the effect of the short-distance repulsion is cancelled by the effect of the long-distance 
attraction). In this case, we no longer need to impose an upper bound on R: it is easy to see that the above 
derivation holds for arbitrarily large R. Hence the result (4.22) is valid for i? —)■ oo in this special case. It is 
easy to see that, for this value of 6, the expression (4.22) is positive, so this special trajectory experiences a 


time delay. Hence, in this special case, we have a gauge-invariant definition of the time delay for a graviton 
propagating in from infinity and returning to infinity. 


4.4 Time advance in EGB: numerical results 

The above perturbative calculation demonstrates that a time advance is possible for tensor-polarized gravi¬ 
tons propagating in the geometry of a small black hole for d > 7. However, several questions remain. As 
discussed above, small black holes are unphysical for d = 5, 6. So for d = 5,6 we will have to study black 
holes which are not small in order to demonstrate that a time advance is possible. Furthermore, we would 
like to determine (for any d) how large the time advance can be. The perturbative result indicates that the 
time advance increases as the impact parameter decreases so we would like to consider 6 as small as possible. 
The lower bound on the impact parameter is 6 ~ L but the above calculation assumes 6 ^ L. Hence to 
determine the largest possible time advance we will need to use a different method. Finally, we discussed 
above the case of special trajectories which experience no net deflection, for which the time delay is finite as 
R —> oo. We saw that, the resulting time delay is always positive when the perturbative calculation is valid. 
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d-8. r^i-0.03^fa 


d-8, rn 




r/^fa b/^Ta 


Figure 2: Effective potential (left) and deflection angle (right) for tensor polarized gravitons scattered by a 
small black hole in d = 8. We set th = O-OSyW, which gives /r « 1.4 x and L « 0.2-yQ'. 


But what about trajectories with b ^ LI Could these exhibit zero net deflection? If so, can they exhibit a 
time advance? 

To address the above questions, we will resort to numerical integration. We will compute numerically 


both the deflection angle (4.111 and the time delay (4.121 for the tensor modes as functions of the impact 


parameter and plot the results for different parameters. In practice, we do this calculation by nsing i?o, the 


minimum value of r, to label the trajectory and determine h from Rq using (4.101. 


We start by calculating the deflection angle for tensor-polarized gravitons by a small black hole with 
d > 7, in the limit of infinite cavity radins. We will compnte the deflection angle A(/)oo as a function of 
the impact parameter b. For d > 8, complicated behaviour arises at small b because the tensor effective 
potential for small black holes has a complicated form with a local minimum and two local maxima - see 
Figure The local minimum corresponds to stable circular graviton orbits around the black hole. The 
maxima correspond to nnstable circnlar graviton orbits. Hence there are three “tensor spheres”. Only the 
local maxima are relevant for scattering of (high frequency) gravitons. As noted above, the deflection angle 
must diverge at impact parameters corresponding to these maxima and this can be seen in in Fig. This 
figure also shows that A(f>ao > tt when Rq lies between the two maxima. 

Figure 1^ shows plots of A(j>oo for d = 7,8 outside the outer tensor sphere. At large b, our perturbative 
results show that > tt although this is not apparent from the plots because A(j)ao — tt is very small. As 

b decreases, our perturbative result shows that becomes negative, as seen in the plots. The plots 

show that A^oo ~ tt decreases to a negative minimnm and then increases, becoming positive as b is decreased 
fnrther. This lies ontside the validity of the pertnrbative calcnlation. Note that there are two valnes of b for 
which A(j)ao = TT, i.e., for which there is no net deflection. The larger of these, with 6 ^ L is encompassed by 
onr perturbative approximation, and we showed above that it gives a trajectory with a positive time delay 
in the infinite cavity limit. However, we will show that the smaller valne of 5 ^ L can give a time advance 
in the infinite cavity limit. 

Figure 1^ plots the time delay for different values of the cavity radius R with d = 7, At large &, the 
time delay is small and positive, but becomes negative as & is decreased, i.e. there is a time advance as 
predicted by our perturbative calculation. The size of the time advance increases as 6 is decreased further 
but at small enough b, the time delay becomes positive again. As expected, increasing R tends to increase 
the time delay. The only trajectories for which this does not happen correspond to the two special values of 
b for which the trajectory does not nndergo a net deflection. Hence as i? is increased, the minimum in these 
plots, corresponding to the largest time advance, becomes more and more localized around the smaller valne 
of b for which A(/)oo = tt. 


^^For d = 8 we only show results for Rq outside the outer tensor sphere. The time delay is positive for Rq between the two 
maxima of the effective potential. 
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Figure 3: Deflection angle for small black holes in d = 7,8. We set rn = O.OSyW which gives fi « 4.5x 10“ 

L « 0.28-y/a in d = 7, and /r « 1.4 x L « 0.2yW in d = 8. The dashed line represents the 

perturbative approximation (4.151. 



Figure 4: Time delay for small black holes in d = 7,8. We set rn = 0.03-ya as before. The solid, dashed, 
dot-dashed and dotted lines correspond to i? = 2.5y/a, R = 5y/a, R = lO-^a and R = 50y/a respectively. 
The black dashed line corresponds to the perturbative approximation (4.22). 


d=7 



Figure 5: Time delay for fixed b (such that A(j>ao = tt) expressed as a function of the radius of the cavity 
for a small black hole in d = 7. We set th = 0.03-y/a as before. The dashed line corresponds to the limit 
i? —>■ oo. 
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d=7 



Figure 6: (Left) Log-log plot of the absolute value of the time advance for an nndeflected geodesic in d = 7 
in the limit i? —7 oo against the mass parameter /i (units a = 1). The dashed line is given by a/i'^ with 
a « 0.09 and k « 0.165. (Right) Plot of k against d. The dashed line is 


Figure shows how the maximum time advance, corresponding to the minimum in Fig. behaves as R 
is increased. As i? —)■ oo we see that the (negative) time delay converges to a finite limit Dao, as expected 
for a trajectory with zero net deflection. We wonld like to nnderstand what scale determines the amplitude 
of Dao- The obvious guess is the scale L and this tnrns ont to be correct. In Fignrej^we plot \Dao\ against 
the mass parameter /i for small black holes in d = 7. From the plot we dednce that the relation shonld be a 
power law: |Iloo| ~ (iii units a = 1). By estimating the valne of k in different dimensions (Fignre|^ we 
obtain numerically n « Recalling that L ^ (since a = 1), we have fonnd: 

\Doo\^L. (4.25) 


One can explain this result analytically as follows. Denote by 6* the value of b close to the tensor sphere 
(i.e. 5* ~ L) for which A^oo = From Fig. [^it is clear that the most of the time advance arises from the 
region r ~ L. For a small black hole, we have / « 1 in this region (see end of section 3.41. This snggests 
that we can calculate D^o for a small black hole by approximating / = 1 in the integral for the time delay. 
If we do this then the integrand now varies only over the scale L and so all quantities in the problem are of 
order L. Hence, by dimensional analysis, mnst be proportional to L. To get an idea of the error made 
in setting / = 1 we note that this approximation eliminates the nsnal GR time delay effect. We can estimate 
the error made by our approximation by estimating the size of this delay as ~ ja. For a 

small black hole L ^ ^Jol hence the error is parametrically smaller than the scale L and therefore negligible. 

In snmmary, we have shown that, for a small black hole with d > 7, there is a tensor-polarized graviton 
trajectory with impact parameter 6 ~ L that experiences no net deflection and, in the infinite cavity limit, 
experiences a finite time advance of order L. 


Finally we consider the cases of d = 5, 6 for which we cannot consider arbitrarily small black holes becanse 
of the failure of hyperbolicityj^ We want to show that a negative time delay is possible for d = 5,6. To do 
this, consider the case in which R = i?o(l + e) with e <C 1. Under the change of variable r = i?o(l + x) the 
integral for At becomes: 


At = 2R, 


da; ^/(i?o(l-kx))^!^ 


/{Rojl -k x))ca{Ro(1 + x))5^ 
i?g(l -kx)2 


-1 


(4.26) 


^®In d = 5 the theory is hyperbolic in the exterior of the black hole for rjjl^/a > \/l + \/2 Ri 1.6, while in d = 6 this happens 
for r£f/\/a > ( W5 ^5 + 2\/6^ — 1 j Ri 0.4 


18 










0.6 

0.4 

Lj 0.2 

Q 

0.0 

- 0.2 

Figure 7: Time delay for tensor-polarized gravitons in d = 5, 6 dimensions. In d = 5 we set = 2y^, which 
gives ^ = 4.5a, L « \.by/a, and we choose R = i^/a. In d = 6 we set rn = which gives = 1.5a^/^, 
L « 1.1^/a, and we choose R = 2y/a. 


d=5, rH=yfa 


d=6. rH=yo 



h/yfa 



We can now expand in powers of a: <C 1 and integrate: 

4i?o-\/e 


At = 


/ 2*/' D t /o (^o) i^o) 

Similarly one can compute the deflection angle: 

__ 4 a / CA{Ro)\/e _ 

fTT~D~\ /o RofiRii) ^oc^(i?o) 

VJ{^o)\J2- FWrTT 

It follows that the time delay is given by: 






(4.27) 


(4.28) 


D = 


4i?o-/e 


i^OC^(Ko) 

caC-Ro) 


(l-v/^/(i^) +0(e3/2). 


(4.29) 


Since cs,c\/ < 1 we see that in this setting the time delay is always positive for scalar and vector modes, as 
expected. For tensor modes this can be negative. In particnlar, in d = 5 we have ct > 1 everywhere and 
thns, for a black hole of arbitrary size we have a negative time delay when R = Rq{ 1 + e). In d = 6, for 
r > L we also have ct > 1- Motivated by this, we compute numerically the time delay in d = 5,6 for valnes 
of R comparable to Rq. The nnmerics confirm that it is possible to obtain a time advance when d = 5, 6 
(Figure]^. We have also studied the deflection angle, which we And is always greater than tt, so zero net 
deflection trajectories do not occur for d = 5, 6. 


5 Time machines 

In this section we will discuss the suggestion that one can exploit the negative Shapiro time delay to construct 
a causality violating spacetime, i.e., a “time machine”, in Einstein-Ganss-Bonnet theories jH]. This argument 
is closely related to arguments applying to any Lorentz covariant held theory with snperluminal propagation, 
some of which have appeared in Ref. HU- which considers various flat space field theories with superluminal 
propagation. Snch time machine constructions have been criticized by Geroch |15| (see also Ref. |16|). In 
this section we will discnss how these criticisms apply to the constructions of Refs Emu. 

Gonsider first the case of General Relativity. There are many solntions of the Einstein eqnation which 
exhibit cansality violation e.g. Minkowski spacetime with a periodically identified time direction. We do not 
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reject GR as a physical theory because it admits such solutions. This is because one cannot “make” these 
time machines starting from initial data. Stated mathematically: such causality-violating solutions are not 
the Cauchy development of any initial dataj^ 

Let’s consider now the argument of Ref. m, which discussed several flat spacetime field theories with 
superluminal propagation. The simplest example is a scalar held with action 


S = - 


1 

2 


d^x 




A4 




(5.1) 


where A is a mass scale and C 3 a dimensionless constant. For this theory, the equation of motion is: 


E{Tr, dn, d^n) 



2 c 3 

A4 


dir ■ dn 


d^d/^n - '^d^nd'^nd^d^n = 0 . 


(5.2) 


For this equation, a hypersurface is characteristic if, and only if, it is null w.r.t. 
metric”: 


G^''{dn) 


1 - 


A4 


the (inverse) “effective 


(5.3) 


Assume that there exists a fiducial inertial frame for which | 9 ^ 7 r| A^j^ Then has Lorentzian signature, 
which implies that the equation of motion is hyperbolic (a frame-independent statement). It is G^'^ that 
determines causal properties of this equation. Following the terminology of the Introduction, we say that a 
covector ^ is timelike iff it is timelike w.r.t. G^" etc. 

Inverting G^'^ gives the effective metric. In the fiducial inertial frame we have 


G^u{dn) 




4c3 „ „ 

-j^d^nd^n. 


(5.4) 


Contracting with a vector gives 

From this it can be seen that the null cones of and 77 ^^ are nested, with the null cone of G^,y inside that 
of 77 ^x 1 / when C 3 > 0 and outside that of rj^i, when C 3 < 0 , i.e., the theory has superluminal propagation when 

C3 < 0 [H]. 

We can now discuss the initial value problem. Given some inertial frame x^, we would like to specify 
initial data ( 7 r, 9 o 7 r) on S = {x° = 0}. Of course S is spacelike w.r.t. rj^'' but for a well-posed problem it is 
necessary that S also be spacelike w.r.t. G^". For C 3 > 0 this is automatic. If C 3 < 0 then this appears to 
restrict our freedom to choose the initial data for 9 o 7 r. But this is not a new restriction: we already imposed 
a restriction on the initial data, i.e., the existence of the fiducial inertial frame. In the fiducial frame, the 
surface x° = 0 is obviously spacelike w.r.t. G^'^. 

If, in some inertial frame, the surface x° = 0 is not spacelike w.r.t. G^" then the initial value problem 
will not be well-posed. In this case, for generic initial data, one would not expect a solution of the equation 
of motion to exist, even locally near x° = 0. One might be able to And a solution for very special initial 
data e.g. if the data is analytic then a solution will exist locally by the Gauchy-Kowalevskaya theorem. But 
this is infinitely fine-tuned: if one perturbs the initial data in a compact region then the resulting data will 

4® More generally, the class of spacetimes of interest in GR is the class of spacetimes that arises as the maximal Cauchy 
development of suitable initial data, where “suitable” depends on the physical situation e.g. one would usually require the initial 
data to be geodesically complete and impose some asymptotic boundary condition e.g. asymptotic flatness. By definition, 
the maximal development is globally hyperbolic and so it can never violate causality. But the maximal development might be 
extendible beyond a Cauchy horizon into a causality violating region, which would capture the notion of formation of a time 
machine. The strong cosmic censorship conjecture asserts that, for suitable initial data, the maximal development is generically 
inextendible. Hence, if correct, strong cosmic censorship excludes time machines because either they can’t be formed or they 
are infinitely fine-tuned (non-generic). 

^®This is probably required for the validity of effective field theory. 
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be non-analytic and no solution can be expected to exist. More generally, the solution does not depend 
continuously on the initial data. 

Ref. [14] argues heuristically that it is possible to construct a time machine when C 3 < 0 by considering 
two lumps of scalar field, well-separated in the direction, which are highly boosted w.r.t. to each other 
in the cc^-direction. So consider initial data at = 0 consisting of two such lumps with a large relative 
boost The problem is that such an initial data surface is not everywhere spacelike w.r.t. m so this 
initial value problem is not well-posed. In general one would not expect any solution of the equations of 
motion to exist for such initial data. So one cannot build a time machine this way. 

One might argue that it is obvious that a time machine could never result from Cauchy evolution of 
initial data since Cauchy evolution will break down when one is on the threshold of forming a time machine. 
So the question we should really ask is whether such a breakdown can occur starting from “good” initial 
data. In more physical terms: if one wishes to employ a large relative boost to build a time machine one 
must specify how this large relative boost will arise from “good” initial data m- 

For the scalar field theory above, “good” means that the initial data surface a:° = 0 should be spacelike 
w.r.t. . Cauchy evolution remains well-posed as long as surfaces of constant remain spacelike w.r.t. 
G^" so a necessary condition for formation of a time machine would be existence of a time T > 0 at which 
the solution remains smooth but the surface x^ = T becomes null w.r.t. G^‘^ at one more more points. This 
would correspond to the threshold of formation of the time machine. 

Can this happen? It is well-known that such behaviour does not occur starting from small initial data, 
i.e., data such that tt and its first few derivatives are small. The solution arising from small initial data simply 
disperses in a similar way to a solution of the linear wave equation So superluminal propagation 

leads to no pathologies in the behaviour of solutions arising from small initial data. Not much is known 
about the global behaviour of solutions of nonlinear wave equations for large initial data. For most nonlinear 
equations, global regularity of solutions is not expected. Solutions can suffer shock formation, i.e., blow-up 
of the field tt (or a derivative of tt) at some time T > 0. See for example Ref. [23] (albeit not for a Lorentz 
covariant equation). As far as we know, it is not excluded that the equation discussed above could have 
large data solutions that evolve to the threshold of formation of a time machine. But, as we have discussed, 
there is no compelling reason to believe that this is the case. 

We now turn to the proposal of Ref. [^ that it is possible to construct a time machine in EGB theory. This 
is done by exploiting the negative Shapiro time delay experienced by gravitons. The proposed time machine 
arises from two high energy gravitons, moving in opposite directions with non-zero impact parameter. Each 
graviton is described by an Aichelburg-Sexl “shock-wave” solutions nnuni. It is assumed that the spacetime 
resulting from the collision is well-approximated by two outgoing Aicheburg-Sexl shock waves. Under these 
assumptions one can argue that there exist closed causal curves in the spacetime. 

One problem with this construction is the use of Aichelburg-Sexl solutions. The curvature of an AS 
solution is a delta-function localized on a null hypersurface (with the amplitude of the delta-function diverging 
on a null line within this hypersurface: this is viewed as the worldline of the graviton). Owing to special 
symmetries of its curvature tensor, this is an exact solution of Einstein-Gauss-Bonnet theory. Now clearly 
one can superpose two such solutions moving in opposite directions (since the spacetime is fiat between 
them) to obtain a solution valid until the two null hypersurfaces intersect. But when these hypersurfaces 
intersect, it is far from clear that there is any sense in which the equation of motion can be satisfied. This is 
because the equation of motion involves products of curvature tensors. Hence along the line of intersection 
of the hypersurfaces, there will be a product of delta functions that cannot be balanced. Therefore it seems 
unlikely that the spacetime can be extended to the future of this intersection]^ 

This problem arises from the fact that an AS solution is singular. So maybe we can solve the problem by 
smoothing out the singularity. As discussed in the Introduction, an AS solution can be obtained by taking 
a limit in which one boosts a black hole solution and takes the boost to infinity whilst scaling the black hole 
mass to zero, keeping the total energy fixed. This suggests that we should consider initial data consisting 
of two small (compared to the GB scale) black holes, moving with high relative boost in opposite directions 

Note that such initial data will not satisfy the condition |9^7r| <C everywhere. 

^^Note that a theory of interacting impulsive (i.e. delta-function curvature) gravitational waves does exist for GR m- 
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with large impact parameter. It would be a difficult matter to construct such data explicitly, solving the 
constraint equations, but there is no reason to doubt that this can be done. 

Now the question is whether this is “good” initial data. As discussed above, black holes with arbitrarily 
small mass don’t exist for d = 5. When d = 6 , small black holes are unphysical because the equation of 
motion is not hyperbolic. So consider d > 7. In section |3.3| we showed that there is a speed limit for small 
black holes arising from the condition that the initial data surface be spacelike. Hence we cannot start 
from initial data describing two black holes with a very large relative boost: such initial data will not be 
everywhere spacelike and hence this data cannot be evolved (or is infinitely fine-tuned), just as for the scalar 
field example discussed above. 

We can attempt to construct legitimate initial data by requiring that the speed limit is respected. Consider 
two small black holes, each of mass parameter /i, boosted in opposite directions with speed v < Umaxi 
separated in the transverse direction by a distance R. Assume that the distance is sufficiently large for the 
gravitational interaction between the black holes to be negligible, i.e. R ^ L. Consider a tensor-polarized 
graviton propagating between the holes. Our numerical results suggested that, for a black hole at rest, the 
maximum time advance \D\ that a graviton can experience is of order L and is achieved for b ^ L. If we now 
boost the black hole in a direction transverse to the motion of the graviton (in order to avoid issues related 
to the length contraction effect), the time advance gets amplified by 7 = (1 — \D\ ^ 7 L. For the 

argument of to work we need the time advance to ‘compensate’ the time taken by the graviton to travel 
between the two black holes, i.e. \D\ ~ R. However for this to hold we would need 7 L ^ R ^ L, that is 
7 3> 1 , which cannot be achieved because of the restriction v < Umax- 

In summary, we have argued that attempting to build a time machine spacetime in EGB theory using 
the method suggested in Ref. [5] will not work. This is because the initial data required is not everywhere 
spacelike (in the sense defined in the Introduction) so the initial value problem is not well-posed: either no 
solution will exist, or it will be infinitely fine-tuned. 
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A Causal structure of a PDE 

We give here a more rigorous treatment of the causal structure of a hyperbolic system, introducing the 
necessary concepts in PDE theory. We refer to EHIISIIIIIZI for the details. For simplicity we will limit our 
discussion to second order PDFs. Consider a second order linear differential operator L on a d-dimensional 
Lorentzian manifold {A4,g), p G Ai. In local coordinates: 

L = {p)dadj 3 + b'^ {p)dj + c{p) (A.l) 

We then define the following: 

Definition A.l. Let ^ S T*A4. The principal part of the differential operator L is: 

PP{L)=a^^ip)d^dp. (A.2) 

The principal symbol of the operator is defined as the contraction: 

= (A.3) 

More generally, for a system of non-linear PDEs, u ■. AA ^ K": 

F/(p, M, du, d'^u) = 0, 1 < 7 < n (A.4) 

one defines the principal symbol as the principal symbol of the linearised system: 
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Definition A.2. The principal symbol of ( A.4[ ) is: 

ot dFi{p,u,du,d'^u) 

aia^dfuJ) 

Definition A.3. Let ^ S T*M.. The characteristic polynomial Q{p,C) is defined asp^ 

Q{P:0 = detP(p,0. 


Definition A.4. We define the characteristic subset of T*Ai (or normal cone at p) as: 

c; = {^&t;m : Q(p,e) = o}. 

Remark. The characteristic polynomial is a homogeneous polynomial in ^ of degree 2n. 
We can now define characteristic hypersurfaces. 


(A.5) 


(A.6) 


(A.7) 


Definition A.5. Consider a hypersurface E C At with normal 1-form £. S is a characteristic hypersurface 

ifg(p,e) = o,VpeE. 

Consider a hypersurface defined as a level set of some function (j). The surface will be characteristic iff 
Q{p, d(^) = 0. This defines a first order non-linear equation for the function (j) know as the eikonal equation. 

Remark. Such characteristic hypersurfaces are generated by bicharacteristic curves. These are defined as 
the curves (a:(A),^(A)) in the cotangent bundle T*Ai satisfying: 


_ dQ{x,^) • dQ{x,0 

"" dx^ 

with the initial condition Q{p,^) = 0, which is preserved along these curves by construction. 
We now want to define the dual of the set C*. 

Definition A.6. The ray cone at p G At is defined as the set 


(A.8) 


= u 

? 6 C; 


•-pA 


(A.9) 


where 

4,5 = {A G TpM : (A,0 = 0} (A.IO) 

Remark. The projections of the bicharacteristics onto the base manifold are called rays and their tangent 
vectors Y (p) = g 4 are called characteristic tangent vectors. The ray cone is nothing else that 

the ‘light cone’ defined by null rays. 

We now define what we mean by hyperbolic systems: 

Definition A.7. The characteristic polynomial Q(p4) is said to be hyperbolic if there exists ( G T*Ai such 
that any straight line with direction such that it does not intersect the vertex ^ = 0, intersects the cone 
QiPi 0 = 0 ill i^al roots, where m is the order of the PDE and n the rank of the principal symbol. 

A system of PDEs is hyperbolic at p G Al if its characteristic polynomial is hyperbolic. We say that the 
system is strictly hyperbolic if, moreover, all the roots are distinct. 

Denote by C* the closure of the subset of T*M defined by such one-forms C. It can be shown that C* 
consists of two convex, opposite cones C* = C'*’+ U C*’“. These are non-empty and dC* C C*. Similarly, 
denote the dual cones in the tangent space by Cp = Cj U C“. In this case we have 4 C Cp. If it is possible 
to continuously distinguish between the convex cones , we say that the spacetime is time-orientable. We 
can finally define the causal structure: 

In a Lovelock theory (or in GR), this definition needs modification because of the gauge symmetry. See Ref. for details. 
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Definition A.8. A vector X S TpM. is causal if A G Cp. A one-form ^ G T* A4 is causal if ^ & C*. 

It follows, in particular, that a hypersurface S is spacelike if its normal one-form is in the interior of C*, 
Vp G E. Thus in discussing the initial value problem we shall prescribe initial data on a hypersurface which 
is spacelike according to this definition. 

Example 1. Consider the class of solutions to Lovelock theory for which the characteristic polynomial fac¬ 
torises as a product of quadratic factors: 

Q{p,i) = [Gf{p)UbT{G‘^2\p)UdY ^... (A.ll) 

It is easy to see that the normal cone corresponds to the union of the null cones of the inverse effective 
metrics, while the ray cone corresponds to the union of the null cones of the effective metrics. The set G* 
corresponds to the intersection of the subsets of the cotangent space defined by G^{p)^a^b < 0. The set Cp, 
instead, corresponds to the union of the subsets of the tangent space for which G/ab(p)A“A^ < 0. In the 
case in which these cones form a nested set (e.g. for Ricci fiat type N spacetimes) then G* (Cp) corresponds 
to the innermost (outermost) cone in TpAi {TpXi). In this case we see that a spacelike hypersurface is one 
which lies always outside of the outermost cone in the tangent space and whose normal vector always lies in 
the innermost cone in the cotangent space. 


B Perturbative calculations 

We give here more details on the perturbative calculation of the time delay and deflection angle. Recall from 
section 4.3 that we want to compute: 

rR 


nn pn 

A(f> = 2b dr CA{'e)r~'^h{r), At = 2 / dr f{r)~^h{r), 
J Rn J Rq 


where we have introduced: 


h{r) = 1 - 


f{r)cA{r)Y‘‘ 


2 \ - 1/2 


(B.l) 


(B.2) 


We want to calculate the above quantities subject to the assumption that i?o is large compared to the black 
hole size in the following sense 

/i , L 




We will assume that the cavity radius is large: 


Uo 


For the time delay we will need to assume that the cavity radius is not too large: 


R r: 


,d-3 


Rq P 


Rq 

~L 


i-l 


(B.3) 


(B.4) 


(B.5) 


B.l Approximation for h{r-) 


The impact parameter is related to i?o by equation (4.10). Using (B.3) we have: 


= R 


2 

0 





R 


d—3 
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0-^)) 


(B. 6 ) 
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similarly, introducing z = Rq Jr (so 0 < z < 1): 


/(r)cA(r) ^ ^ M l^^d-3 _ 2a/3A ,d-i'^ j 


Rl 


i?i 


d-3 


Rl 


And hence: 


Me HI+ ((!-»«) 

((^ 


2aPA, 


Kq 


and the ellipsis denotes terms of order O 


Tci — 3 ' D<i—1 


B.2 Approximation for the deflection angle 

Changing the integration variable to z: 

b 

A(j) = 2 — / dzcAh 


i?, 


0 jRo/B. 


In our approximation, denoting by ... terms of order O 


,<1-3 ^ r,d-l 


, we obtain: 


[}^Rd-3[2^ (l-z2) l-z2 


Which yields: 


Acj) = 2 arccos ( ^ ^ 


RJ ' R'^-^ 




+ ... 


where we have defined: 


j=f' d^d-z^r'i^ (] 

JRo/R \‘ 


1 2 (1 - z‘^~n a^A fz^ + z‘^+^ - 2 z' 


(l-z2) Rl 


^d—1 


1 -Z 2 


Using the large cavity radius assumption (B.4) we have J = Jq + 0{R^/R^), where 
. 2a/?A(d-2)\ T (f) tt /I a/3^\ 

Jo = H(<1-1) (1 -- 5 (5 -. 


Moreover we have that 


and hence 


0 I -^0 A 0^^ n f 

2 »ccos|- =,- 2 -+ 0 U 5 


^,1 af3A\ fRo\,,r 2aPA{d-2)\ r(^) , 

2 (.-^]arccos[-j+ 2 J = 0 f(d-l) ^ j 21 ^ + ^ 


.2 Rl ) 

We can then conclude that 


i?o 


A<(. = TT - 2—+ 


+ 0 


R R^ 


Rn 


■ / 2aPA{d-2)\ r(|) 


^d -1 


i? ' V<”^ 


(B.7) 

(B. 8 ) 


(B.9) 

(B.IO) 

(B.ll) 

(B.12) 

(B.13) 

(B.14) 

(B.15) 

(B.16) 
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B.3 Approximation for the proper time 

Denoting again by ... terms of order O + ^d-i ^ ^, 


we have: 


= (1 - (i + + 


1 (1 - 2aPA (1 - 


2 1 - z2 


i?2 1-^2 


from which: 


where: 


At = 2jR^-m + 2 


0 ' rid —4 

/to 


/ + ... 


I = 


Jro/R \ 


l(l_^d-3) a/3^ (1 _ 

2 (l-z2) 




For large cavity radius, i.e., (B.4|, we have: 


At = 2i?-^+2^/o + ofi?o 


1 — 


-d-l 


i?0 


R ' Rt^ 


In = ^\l- 


where 


Plugging back in the above we obtain 
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2apA{d-^)\{d-l){d-i) r(f) 
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Rl 

R 


(B.17) 

(B.18) 

(B.19) 

(B.20) 

(B.21) 


(B.22) 


Finally, note that in this approximation we have f{R) = 1 + 0{ j^a_3 ) and ^ = {fJ,/Rl ^){Rq/RY ^ 

which is negligible to the order of approximation used above, hence At « At. 


B.4 Approximation for the time in Minkowski 

We have: 

M f _ Q:/ 3 a \ 

1,1,2 Rl ) 
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We now use the condition (B.5) that the cavity is not too large. This ensures that the last term above is 
small and we obtain: 


AtMink -2R + (V^{d 1) (^1 ^ ^ . 
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B.5 Approximation for the time delay 


Finally, putting the above results together, we have determined the time delay under the conditions (B.3l, 

ra, ra: 


D = 
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d-i 
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/ 2a/3^(d-4) \ (d-1) r(|) 

V Rl (d-i)y (d-4)2^(^^) 


+ 0 




(B.25) 
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